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Abstract 

We study a system of rods on Z 2 , with hard-core exclusion. Each rod has a length between 
2 and N. We show that, when N is sufficiently large, and for suitable fugacity, there are several 
distinct Gibbs states, with orientational long-range order. This is in sharp contrast with the 
case N = 2 (the monomer-dimer model), for which Heilmann and Lieb proved absence of phase 
transition at any fugacity. This is the first example of a pure hard-core system with phases 
displaying orientational order, but not translational order; this is a fundamental characteristic 
feature of liquid crystals. 


1 Introduction and results 

In 1949, Lars Onsager proposed a theory of the isotropic-nematic phase transition in liquid crystals, 
which relied on the following simple heuristics m- Picture each molecule as a (very) long, (very) 
thin rod. There is no energetical interaction between the rods, except for hard-core exclusion. 
Since at low densities the molecules are typically far from each other, the resulting state will be 
an isotropic gas. However, at large densities it might be more favorable for the molecules to align 
spontaneously, since the resulting loss of orientational entropy is by far compensated by the gain 
of translational entropy: indeed, there are many more ways of placing nearly aligned rods than 
randomly oriented ones. 

This is probably the first example of an entropy-driven phase transition. It shows that an 
increase of entropy can sometimes result in an apparently more ordered structure, hence the often 
used expression “order from disorder”. 

In spite of the obvious physical relevance of such issues, rigorous results are still very scarce. 
The only proof of such a phase transition has been given (as a side remark) in jj for the following 
simple model: The rods are one-dimensional unit-length line segments in R 2 , with two possible 
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orientations (say, horizontal and vertical); a configuration of N rods (N is not fixed) is specified by 
a family (x,a) £ R 2Ar x {—1,1}^, where (x 2 k-u% 2 k) is the position of the middle of the £:th rod, 
while (jfc represents its orientation. A configuration in a subset V of R 2 is admissible if all rods 
are inside V and are disjoint; let Cy denote this event. Then the measure describing the process 
is v\(-\ Cy), where v\ is the product of the Poisson point process in R 2 of intensity A > 0 and 
the Bernoulli process of parameter A. The main result is then that, in the thermodynamic limit 
V y R 2 , there are (at least) two limiting Gibbs states with long-range orientational order, for all 
A large enough. This model has a very special feature though, namely two horizontal (respectively, 
vertical) rods have O-probability of intersecting under v\. This is a considerable simplification, and 
therefore this result does not provide any information on the much more interesting case of rods 
of finite width. 

The other rigorous results are concerned with lattice versions of this problem. Heilmann and 
Lieb proved in a classical paper [2] that there is no phase transition in the monomer-dimer model, 
in the sense that the corresponding free energy is always analytic. This model is defined as follows 
(on Z 2 , their result holds more generally). Let V be a finite subset of Z 2 . For x,y £ Z 2 , we 
write x ~ y if \x — y\ =1. Let £\> = {{x,y} C V : x ~ z/}; for e, e' £ £y we write e ~ e' if 
e n e' / 0. The state space is given by f2 = {0, l} £v . A configuration to £ is admissible if 
e ~ e' ==> min (u> e ,u> e i) = 0. The probability of a configuration ui is then given by 

/xaM oc 

is admissible} ^ ? 

where |w| = an< ^ > 0- Informally, when a pair of sites e is such that uj e = 1, then the 

two sites are occupied by a dimer; a configuration is admissible if no site belongs to more than one 
dimer; A is the dimers’ fugacity. 

An alternative approach to this model was then discovered by van den Berg GU Using dis¬ 
agreement percolation methods, he was able to give a very simple proof of the much stronger result 
that this model is in fact completely analytic, in the sense of 0. This paper is also interesting 
in that it clearly points out the very special nature of dimers. Indeed, it would be impossible to 
push the analysis to arbitrary values of fugacities, were it not for a magical property of dimers: A 
monomer-dimer model on a graph G is actually equivalent to a pure hard-core gas on the line-graph 
of G. 

Several other similar models have been introduced (see, e.g., 0EH3), in which existence of 
orientationally ordered states has been proven. All these models, however, share the same defect, 
namely the ordered states also automatically display long-range translational order, i.e. they are 
perturbations of periodic configurations. Thus they really can’t be considered satisfactory models 
of liquid crystals, since a central characteristic of the latter is the liquid-like spatial behavior in 
the ordered states. In order to solve this problem, Heilmann and Lieb [21 proposed five different 
models of hard-core particles (actually dimers). For these models they proved the existence of 
long-range orientational order at low temperatures, and gave quite plausible arguments in favour 
of the absence of long-range translational order. These, however, were not pure hard-core models, 
since an additional attractive interaction favouring alignment of the dimers was introduced, and 
thus the question of whether pure hard-core interaction can give rise to such phases was left open 
(actually, Heilmann and Lieb even stated that it was “doubtful [...] whether hard rods on a cubic 
lattice without any additional interaction do indeed undergo a phase transition ”). 

To the best of our knowledge, these are the only rigorous results pertaining to this problem. It 
would be extremely desirable to prove the existence of a phase transition in the monomer-fc-mer 
model (replacing dimers above by /c-mers, i.e. families of k aligned nearest-neighbor sites), for large 
enough k. This seems rather delicate however, and in this work we concentrate on another variant 
of the monomer-dimer model, with only hard-core exclusion and for which it is actually possible to 
prove existence of phases with orientational long-range order and no translational long-range order; 
actually it is also seems possible to treat the three-dimensional case, which presumably would lead 
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to a counterexample to the above claim. We hope to return to the monomer-femer problem and 
to the case of higher dimensions in the future. 

Our model is defined as follows. We call rod a family of k, k £ N, distinct, aligned, nearest- 
neighbor sites of Z 2 a ferod is a rod of length k, and we refer to 1-rod as vacancies. Let V C Z 2 ; a 
configuration ui of our model inside V is a partition of V into a family of disjoint rods. We write 
Nk(co) the number of k- rods in ui. The probability of the configuration u> is given by 

Pq,N,v{u) OC l{JV fe (w)=0, \/k>N}(2q) Nl( ' UJ ' > q^ k = 2 Nk ^ , (1.1) 

where q > 0 and N £ N. Informally, only rods of length at most N are allowed; the activity of 
each rod of length at least 2 is q, and is independent of the rod’s length; there is an additional 
activity 2 q for vacancies. 

Remark 1.1 The activity of vacancies can he removed at the cost of introducing an additional 
factor ( 2q)~ k for each k-rods (k ^ 2). 

Our main task is the proof of the following theorem, which states that for large enough N, 
there is a phase transition from a unique (necessarily isotropic) Gibbs state at large values of q to 
several Gibbs states with long-range orientational order at small values of q, but no translational 
order. This is thus the first model, where such a behavior can be proved. 

Theorem 1.2 1. For any N ^ 2, there exists qo = qo(N) > 0 such that, for all q ^ go there is 

a unique, isotropic Gibbs state. 

2. For any q > 0 sufficiently small there exist Nq = No(q), such that for all N ^ No there are 
two different extremal Gibbs states with long-range orientational order. More precisely, there 
exists a Gibbs state p q N such that 

(iq, N (0 belongs to a horizontal rod ) > 1/2. (1.2) 

In the sequel we shall refer to the infinite volume Gibbs state p q N as to the horizontal state. 
By symmetry the n/2 rotation of the latter gives the vertical Gibbs state p qN , which would 
statistically favour vertically oriented rods. 

A funny consequence of the techniques we develop in order to prove Theorem ll.2l is the following 
result on a sampling of infinite volume horizontal and vertical states by the shapes of the family 
of finite volume domains: 

Theorem 1.3 Let k = (fe, £ 2 ) be two natural numbers. For n = 1,2,..., consider lattice boxes 

Vjf = [—km ,..., km] x [—fen,..., fen], 

and let p q N y k be the finite volume Gibbs state specified in O- Then if q and N satisfy conditions 
of 2) of Theorem \l.‘A 

Hm k'n N V— = Fq,N k l > k 2 

and 

lim l^a N V— = Vq,N */ k l < k 2 

71 —KX) Hi ly i v rL 'll 

Theorem o is proved by showing that, for N large enough, the model defined above is a 
small perturbation (in a suitable sense) of the “exactly solvable” case N = 00 . For the latter, the 
theorem takes the following form. Let q c = 1/(2 + 2y/2). 

Theorem 1.4 1. Let N = 00 . For all q ^ q c there is a unique, isotropic Gibbs state. 

2. For all q < q c , there are (at least) 2 different extremal Gibbs states with long-range orienta¬ 
tional order. More precisely, there exists a Gibbs state p q such that 

p q { 0 belongs to a horizontal rod ) > 1/2. 
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2 An exactly solvable case: Proof of Theorem 11.41 

In this section, we show that the model obtained by setting N = oo is actually exactly solvable, 
since it can be mapped on the 2D Ising model. 

We suppose that our system is contained inside a square box V of linear size L; we suppose 
that we have periodic boundary conditions. We want to partition V into two disjoint subsets 
corresponding to the regions occupied by horizontal and vertical rods respectively. This can be 
done easily once we have said what we do with vacancies. The trick is to split vacancies into two 
species, horizontal and vertical. Doing so, starting from a configuration ui of our original model, 
we obtain a family of 2 iV| different configurations 2*, i = 1,..., The probability of each 

such configuration is then taken to be 

/*,,*,v(£5) cx l{7Vfc(u>)=0, \/k>N} 2- Vl(S) (2 q) Nl ^ = gESLi W(S). 

We can now partition V = 14 V V v into two disjoint subsets. Once these subsets are fixed, the 
problem is reduced to the study of one-dimensional partition functions; indeed each maximal 
connected horizontal piece of 14 can be filled by horizontal rods independently of what choice is 
made for the rest of the configuration, and similarly for vertical pieces of 14- 

In the N = oo the one-dimensional partition functions could be computed exactly: 

z “ = E ( n :') « ,+1 = i a+a+«>” ■ 

where if 40 is the ID partition function in a box of length n ^ 1. Now observe that the exponentially 
decreasing term (1 + q) n is actually irrelevant, since its total contribution to the weight of a 
partition 14 V V v is (1 + qy^ , and is therefore independent of the partition. We thus see that 
this model possess the remarkable property that all its ID partition functions are actually equal 

to e -4 ^ = q/( 1 + q ), and therefore independent of n. It is then very easy to compute the total 
weight of a partition: 

weight(I4, V v ) oc e _2/3 ^l , 

where 7 = (71 ,...,y m ) is the set of contours of the partition, i.e. the set of all bonds of the 
dual lattice intersecting a bond between two nearest-neighbor sites belonging one to 14 and the 
other to V v ; I7I is the total length of the contours, where the two components of the partition are 
reinterpreted as the components occupied by +, respectively. — spins. 

One thus observes that the weight of partitions are the same as those of the corresponding 
configuration of the 2D Ising model at inverse temperature /?, in the box V with periodic boundary 
conditions. Now notice that /3(q c ) = ^ log(l + y/2) = /3 C , the critical inverse temperature of the 
2D Ising model. It immediately follows that for q ^ g c , the corresponding Ising model is in the 
high-temperature phase, and therefore possesses a unique Gibbs state. Statement 1 of Theorem ll.4l 
follows immediately from the symmetry of the latter Gibbs state. 

To prove statement 2. requires only a simple additional argument. For a given collection of 
rods 2, let us denote by Z ^ the number of sites of 14 containing vacancies, and = |I4| — Z h; 
similarly introduce N v and Z v . When q < q c , the Ising model is in the low-temperature region; 
consequently, 

E, 9 , v [||I 4MK||] = Eia ing „3,y[|y^ CT a |] > cL 2 , 

xev 

with c > 0. Since fi q y(0 belongs to a horizontal rod ) = L~ 2 v [A4] , the conclusion now 
follows easily from 

||I4|-|K|| ^ |V h - Ny\ +\Z h ~Zy\, 
and v [|i4 — Z v |] < CL, by the Central Limit Theorem. 
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Remark 2.1 A lot of additional information (e.g., on the critical behavior) can be extracted from 
this mapping to the 2D Ising model. We refrain from doing that here, since this is quite straight¬ 
forward. .. 

Remark 2.2 As it was pointed to one of us by Lincoln Chayes, in N = oo case the techniques of 
reflection positivity enable to treat a more general situation when the rod weights are given by 

\Nx{u) qN k (ui). 
k 

In the above notation the case we consider here corresponds to a specific choice A = 2 q. However, 
the reflection positivity argument does not go through when there is a finite collection of admissible 
rod lengths, N < oo. 

3 Asymptotics of one-dimensional partition functions 

Our next step is to show that the model with finite (but large) N is actually a small perturbation 
of the exactly solvable model analyzed in Section [3 The idea, which is described in details in 
Subsectionis to replace all the ID partition functions by their limiting values (for n —>■ oo), and 
to expand the error term. To be able to control this expansion, we need a very good control on the 
speed of convergence of these ID partition functions. This is the aim of the current subsection. 

3.1 The setup. 

We shall consider here a general case of non-negative rod activities {fk} which we shall view as a 
perturbation of the geometric distribution, 

fk = QP 1 *- 1 +e fc ; k = 1,2,..., (3.3) 

where p + q = 1 and the activities {fk} are normalized to furnish a probability distribution, that 
is 

E e * = 0 ( 3 - 4 ) 

k 

The important assumptions are those on the smallness of the perturbation sequence {efc} with 
respect to the background geometric distribution {qp k ~ 1 }'. 

Assumption A1 There exist <5 < oo and p £ (1 ,p~ 1 ] such that 

\e k \^Sp~ k , A: = 1,2,.... (3.5) 

Assumption A2 There exists a > 0 sufficiently small such that, 

<5 < a(p — l) 2 . 

Assumption A1 is an essential one. On the other hand, Assumption A2 is more technical and it 
merely reflects an intended compromise between giving a relatively simple proof and yet generating 
a whole family of examples where the entropy driven phase transition takes place. Notice that since 
p < 1/p, assumption A2 in fact implies a bound on 6 in terms of q: 


Given {fk} as in (13.311 above we use it to set up the renewal relation: 

n 

go = l and g n = ^ fkg n -k n = 2,3, — (3.6) 

fe=l 
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Define the generating function of the {/&} sequence as 

OO £ 

m) = E f^ k = 7 ^ 7 +E e ^ k = Q(o+ E ( 0 - 

fc=l ^ fc 

Above Q is the generating function of the geometric distribution { qp k ~ 1 } and, accordingly, E is 
the generating function of {e*,}. 

In the sequel we use the notation 


D r .(x) = {z £ C : \z — x\ < r} 


for an open complex ball of radius r centered at x. By Al, F is analytic on 1D> P (0). 

The generating function G of the {g n } sequence is defined and analytic in {z £ C : |z| < 1}. 
By the usual renewal theory, 


lim g n 

n —>oo 


l 


1 A 

i/^ + Efc^fe 


(3.7) 


3.2 The representation formula 

For every v £ (0,1), 

= _L I f dC _ d£ ) 

9n 9 2ni f Du{0) \e(i-F(o) 

(3.8) 

J_/ d£ f (F(g) - 1)-¥'(!)(£-1) ] 

2nif D A n \ (F(0-1)K-1)F'(1) J 


Now, 


As a result, 


(F(0 - 1) - F'(l) = {(Q(0 - 1) ^ Q'(l)} + {E(0 - E'(l)(£ - 1)1 


27rzF'(l) (g n - g) 


I dg f (Q(£) — l) — Q^i)^ — 1) 

JdA h \ (F(0-!)(£-!) 


r d^ f E(Q-r(i)(g-i) l 
JdA u l (F(0 — i)(£ — i) J 


+ / D Ju.<f ) - 

(3.9) 

On the other hand, since the geometric distribution {qp k ~ 1 } generates (via the renewal relation) 
a constant sequence {<?}, 


a r A f (Q(Q-l)-Q'(l)(£-l) 

JD,A n \ (F(0-1)«-1) 


/ dg r (Q(o-i)-q>m-i) i _ n 

JdA h l (Q(^) - !)(€-!) J" 


Subtracting the above expression for zero from the first term on the right hand side of (EHil) we 
obtain 



(F( 0 - 1 )(Q (0 -!)(£-!) 


A 


D 


, F U2<0 ' 


(3.10) 


Thus, with Ui and U 2 being defined as in (TOll and 633 above, the representation formula 
reads as 


27tzF'(1) ( g n 


g) = £ p{Ur(0+U 2 (0}. 
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Assume that it so happens that both Ui and U 2 are analytic in an open neighbourhood of 
]D>ij(0) for some R > 1. Then (Id.811 implies that 

| 9 „- 9 |< (max|Pi(0| + max|U 2 (0|) (3.11) 

We shall represent Ui and U 2 as ratios of two analytic functions, In Subsection Id.Ml we derive a 
lower bounds on the denominators, whereas in Subsection ei we derive the corresponding upper 
bound for the numerators. Eventually we shall pick R = (1 + p)/2 and the target bound on 
(max|£| =fl |U(,j)| + max| 5 | =fl |U 2 (£)|) is formulated in Subsection Id.51 Finally the case of uniform 
rod weights is worked out in detail in Subsection ED 


3.3 Lower bounds on the denominators 


By A1 the function 


m -1 
£-1 


1 

1 


E(0 a 

£-1 


1 

1 


+ V(£) 


is analytic in B p ( 0 ). 

Let us pick 77 G (0,p— 1), later on we shall settle down with the choice 77 = (p — l)/2, but in 
principle all the estimates below could be further optimized. Since 1 + p < p < 1/p, 


inf 

£eB 1+r ,(o) 


1 

1 




1 1 

1 + (1 + r])p ^ 2 ' 


(3.12) 


It, therefore, remains to derive an appropriate upper bound on |E(£)/(£ — 1)| = |V(£)|. There are 
two cases to be considered: 


CASE 1. £ € Di + 7 ) (0) \ 1^(1). Then, by Al, 


|V(£)| < - 



(5(1 + 77 ) 

7 l{p~ (1 + v))' 


(3.13) 


CASE 2. £ G ID^(l). Since E(-) is analytic in B^(l), 

00 00 k / 7 \ oo 

E (£)= 5Z efe ^ fc = 5Z efe 5Z ( 7 ) (£ _ 1 ) ! = 5Z £; ^ _ 1 ) 1 ’ 

1 1 7=0 ' ' 1 

where we have used e*, = 0 and, accordingly, have defined 


<5 



In view of the assumption Al, 


6 i <«!>-*() 


k=l 


= 5 


Sp 


Consequently, 




(1-1 /p) 1+1 (p- i) 

s P 


(p- i)"*- 


l—l 


(P- 1) V (/»-!)' (P — 1)(P — (1 + r ?)) 


(3.14) 


(3.15) 


whenever £ gD,(l). 

Pick 77 = {p— l)/2. Then the right hand sides of both 13.1311 and 13.151) are bounded above by 
25(1 + p)/(p— l) 2 . Only at this stage we evoke assumption A2: under an appropriate choice of a 
the latter expression is as small as desired, say less than 1/6. In view of 13.121) . 13.1 dll and 13.1511 
we, therefore, conclude: 
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Lemma 3.1 Assume A1 and A2, Then, 


min 

feO(i+ P )/ 2 (o) 


no -1 

c-i 



Finally, by direct computation: 


Q(0 -1 


1 

C-i 


1 -p€ 




1 1 

1 + (1 + rj)p ^ 2 ' 


(3.16) 


(3.17) 


3.4 Upper bound on the numerators 

We continue to employ the notation of the preceeding subsection. In particular, 

V(0 = ^ an d ^-E'(1)=V(0-V(1). 

Since by 13.41 , e k = E(l) = 0, ¥ is analytic on D p (0). As in Subsection 13.31 pick 77 £ (0, p — 1) 
and consider the following two cases: 

CASE 1. £ £ Oi +r) (0) \ 10^(1). By (I7TTTT1) and )) 

V(0~V(1) Sp <5(1 + 7?) 2 Sp 

C-l ^ v(p- !)(p- 0 + v)) V 2 (p- 0 + 0) " V 2 (p~ 0 + 0)' 


CASE 2. £ e 15^(1). Since V is analytic on D^(l) we , employing the notation of Subsection 13.31 
estimate: 


v(0-v(i) 


€-1 



1 =2 


Sp 

(p-l) 2 (p-(l+77))’ 


where we have performed a straightforward series summation bounding |e)| as in i rm . 
Picking 77 = (p — 1 )/2 and R = (1 + p)/2 we infer: 


Lemma 3.2 Assume Al, Then, 


max 

€£D h 


E(Q-(^-i)r(i) 

(£-l ) 2 


2 Sp 


(3.18) 


On the other hand, 


1 / Q( 0 -i 
£-1V C-1 



p 

9(1 -pO 


Since by assumption Al, maxi^i = jj|E(^)| ^ Sp/(p — 1), we arrive to the following bound for the 
numerator of U 2 : 


max 

\S\=R 


(Q(C)-i-Q'(i)(£-i))E(e) 

(e-i ) 2 


P _ Sp 

q{ 1 - Rp) (p - 1 ) 


2 Sp 

^ p(p- i) 3 ’ 


(3.19) 


3.5 The target bound on \r n \ = \g n — g\ 

As before set R = (1 + p)/ 2. By the estimates of Lemma till and Lemma f3 .21 


max |Ui(£)| + max |U 2 (0I < 
I? s: R\ \i^R\ 


66p(2+p) 

p{p~ l ) 3 ‘ 







































Finally, F'(l) = 1 /q + J2 k £ fc ■ By assumption Al, 


^ i r p - i q 

By the scaling relation between p and S (assumption A2) the right hand side above 
particular, 

9 = 9(1 + °(l))j 

as it now follows from E3- 

Substituting the above estimates into lETTTll : 

Theorem 3.3 Assume Al and A2. Set R= (1 + p)/2. Then, 

1 1 , , ^ 12qSp(2 + p) 12q5p(2 + p) 

\r n | = \9n ~ g | ^ -7-7T75— R = 


P(P~ l ) 3 


p(p-l ) 3 V! + P 


In particular for every v < (p — l)/4, 


5>ia + o" - Eb - fc id + o" < 


48p(2+p)(l +p) a , 

€ - p(p-l ) 4 - 69 = C2 ^P>° 9 ' 


is o(q). In 
(3.20) 


(3.21) 


(3.22) 


3.6 Uniform rod activities 


Let the rod activities be given by 


fk 


q-,k = l,...,N 
0 ; otherwise 


(3.23) 


Above q = (Y^i P k 3 ) 1 = (1 — p)/( 1 — P N )- Thus, q — q = qp N /(1 — p N ). In other words the 
sequence of weights {/*,} in 13.2311 corresponds, in the notation of 13.311 . to 




qpN 1 *- 1 :k<N 


P 


1 -p N 

- gp k ~ x 


;k>N 


(3.24) 


Without loss of generality we may assume that q < 1/2. Then for each N fixed the weights {e^} 
satisfy assumption Al with , for example, 

P= + “) / 2 = 1 + 2(1-q) and 6 = 5n ^ = i 1 ~ f) ■ ( 3 - 25 ) 

Of course, assumption A2 will be also satisfied for such choice of p and <5 as soon as 

M?) = (1 - §)" « o (p -if = * (jjjUy) ■ 

For the value of p related to q as in (TT25ll set 


ci(q) = ei(q, p) and c 2 =c 2 (q,p), 

where ci and c 2 are the universal constant which appear on the right hand sides of and 

E23. Let us reformulate the claim of Theorem 13.31 as applied to the case of uniform rod activities 
(with the scaling choice 13.251) in mind): 
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Lemma 3.4 Let q < 1/2 be fixed. Then there exists No = No(q) such that for every N ^ No, the 
uniform rod weights {fk} in 113.2311 generate the renewal sequence {g n } which satisfies: 

\r n \ = |ffn ~ g\ < ci(q)S N (q) ^1 + ^ , (3.26) 

where g was defined in <E3. Moreover, for v = (p — l)/4 = < 7 / 8(1 — q), 

5^W(l + *') n ^ c 2 (q)6 N (q)g. (3.27) 

n 

4 Perturbation theory 

In this section V is the lattice torus of a fixed (large) linear size L; V = Z 2 /mod (L). Notice, 
however, that all the estimates below do not depend on L. 

4.1 Super-contours 

We proceed similarly as in the proof of Theorem II .41 We first split vacancies into two families, 
and partition the box V into the two disjoint sub-boxes 14 and V v containing the horizontal, resp. 
vertical, sites. Associated to this partition, there is a family of one-dimensional boxes A = (A;), 
each of which is either a horizontal “segment” in 14, or a vertical “segment” in V v . The weight 
of the partition can then be expressed as a product over all A G A of the corresponding one¬ 
dimensional partition functions 3 ja| — ■ Contrarily to what happens in the case considered in 

Theorem Q these one-dimensional partition function do generally depend on the length of the 
corresponding box A. However, as we have seen in Section |3 these partition functions approach 
their limiting value g rather quickly, provided we choose N large enough. It is therefore convenient 
to expand them around this limiting value: 

AeA AeA v y 7 

We want to use this expansion in order to obtain a perturbation of the pure Ising model which 
appeared in the case of Theorem rrn Let us denote by 7 = ( 7 *) the family of Ising contours 
appearing when interpreting 14, resp. 14, as the region occupied by +, resp. —, spins. We can 
then associate to each of these contours a weight 10 ( 7 ) = e~ 2 ^\ where we have set e~ 2f} = yd?; 
this allows us to write simply 

n.,=n e-^ 1 ■ 

AEA 7^7 

We would like to encode all the information from the partition into these contours; in order to 
do this, we suppose that these contours come with a “color”, i.e. each contour 7 carries the 
information on which of the two sets 14 or W belong to which side of the contour. Of course, there 
is then a compatibility condition on these contours (in addition to their being disjoint): the colors 
must match. 

We also introduce the set of excited intervals / = (J,) C A, and associate to such objects the 
weight w(I) = (<771 — g)/g. Using this we can write 

n( i +2 ^) = En »<«). 

AeA v » / /ca 1 el 

Of course, since our colored contours 7 contain all the information on the partition, the family A 
is actually completely determined by the contours. 
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Figure 1. A configuration of super-contours on the torus V. There are four super-contours. 
Notice that T if of horizontal type, but also possess some interior components of horizontal type 
(shaded in the picture). There are two super-contours winding around the torus; the contribution 
of configurations containing such super-contours being negligible when the box is large, as shown 
in Subsection rm they can actually be neglected. 


We now introduce our basic notion of super-contours, which are maximal connected components 
of (colored) contours and excited intervals (saying that an interval is connected to a contour if at 
least one of the extremities of the interval belongs to the contour). We denote the family of super¬ 
contours by £ = (]?j). The weight u’(T) of a super-contour T is then naturally given by the product 
of the weights of the contours and excited intervals it encompasses. We therefore finally obtain the 
following expression for the weight of the total partition function of our model: 

z q ,N, L = e n w ( r ) * 

r rer 

where the sum is taken over all compatible families of super-contours, i.e. those resulting from a 
partition V = 14 V V v in the way just described. 

At this stage, it is not possible to apply a simple Peierls argument in order to control our model. 
Indeed, our super-contours are colored, and even though there is a symmetry in our model (under 
a simultaneous rotation by 7r/2 and exchange of horizontal and vertical sites. On the other hand, 
there is also a fundamental asymmetry: The shape of a region generally strongly favours one of 
the two species. This forces us to use the general strategy of the Pirogov-Sinai theory, which turns 
out to be quite simple in our case, due to the fact that, because of the above-mentioned symmetry, 
the free energies of the two phases are necessarily equal, and thus we are not required to add a 
suitable external field to reach phase coexistence. 

The basic idea of the Pirogov-Sinai theory is to expand the partition function only over external 
contours, and introduce new weights, which reduce the compatibility condition to something of 
purely geometrical nature. 

However, since we work on the lattice torus V = Z 2 /mod(L), the notion of exteriour of a 
contour is ambiguous. One way to mend the situation would be to fix a distinguished site, say 0, 
and to declare it to be “a point at infinity”. On the other hand all our computations below are 
based on relatively crude combinatorial estimates which take into account local graph geometry 
of Z 2 , but not the global topological structure of V. Consequently, we shall from the start ignore 
(necessarily long) winding super-contours and then simply notice that should we use the “point 
at infinity” definition of exteriour, the analog of below would anyway render long winding 

contours improbable. 
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For any non-winding super-contour T the exteriour of T is defined in a straightforward fashion 
and, accordingly, the type of a non-winding super-contour T will be declared to be horizontal or 
vertical if such is the colour of its exteriour. 

Thus, let (respectively Sf ) be the set of all non-winding horizontal (respectively vertical) 
type super-contours on V = Z 2 /mod (L). Of course, S £ and S'f are related by the 7 t /2 rotation 
symmetry: IfT G S'f, then 9 V / 2 T G Sf. The interiour int(r) of coloured super-contours T G 
is also coloured and in the sequel we shall write 

int(T) = inth(T) U int v (r) 

for the horizontal and vertical parts of int(T). By the 7 r/ 2 -rotation symmetry of V, restricting to 
external contours of horizontal type (“h-type”), i.e. with their exterior colored as horizontal, yields 
exactly one-half of the full partition function: Using the notation 5^' cxt f or se t 0 f a ll collections of 
compatible external contours from S £ and, accordingly, for for set of all collections of compatible 
contours from S } f, we can then write 

Z q ,N,L = 2 E II Z int h (r) Z int v (r) 
reS h, ex t rer 

= 2 E n “(T) ^t(r) 

FeS h, ex t rer 

= 2 E 

res£ rer 

where the new weights are given by 

^ ( r) = w (r)^£). (4.28) 

Z int v (r) 

Notice that in the last expression the sum is over all compatible families of h-type super-contours; 
in particular, the compatibility condition is now purely geometrical. 

4.2 Cluster expansion 

The next step is to show that the new weights are still under control. We first need a bit of 
terminology. Let us denote by T ^ T' the relation ‘T is incompatible with T'”. A cluster is a 
family C of super-contours which cannot be split into two disjoint families C\ and C 2 such that all 
pairs T G Ci and T' G C 2 are compatible. We also write C ^ T if there exists T' G C such that 
T' 7 ^ T. Finally we write |T| for the total length of all the contours and intervals forming T. We 
want to be able to use the following classical sufficient condition for the convergence of the cluster 
expansion 0 : 

Lemma 4.1 Suppose that, for some small a > 0, 

E e 2 “l r 'l|w(r / )| <a|r|, (4.29) 

r' : r'/r 

for each T. Then Z q n,l 7 ^ 0 and there exists a unique function <f> T on the set of clusters such that 

log z qNL = e $T ( C )- 

ccv 


Moreover, 


E |$ T (C)|e a|c| < a|F[. 


( 4 . 30 ) 



We claim that the weights w indeed satisfy gut once q is chosen to be small enough and then 
N is chosen sufficiently large. The argument comprises two steps: First we shall check that gut 
holds for the weights w. Next we shall argue that the conclusion of Lemma f4. II for the weights w 
actually implies the validity of 11X2311 for the target weights w for a possibly smaller value of q and 
larger values of N. 

Lemma 4.2 There exist a > 0, q > 0 and No = No(q) such that 

Y e 2 a|r,| K r ')l ^ a |r| , (4.31) 

r': r'/r 

for every N ^ Nq. 

Lemma 4.3 There exist a > 0, q > 0 and No = No(q) such that (14.2911 holds for q and for every 
N > N 0 . 


4.3 Proof of Lemma 14.21 


A convenient way to over-count 


is as follows: Pick 


Y e 2a i r 'iHr')l 

r' : rvr 


2 a = 


P- 1 


(4.32) 


16 32(1 — q) 

Any excited interval I of a super-contour T connects two dual bonds b and b' which belong to Ising 
contours 7 and, accordingly, 7 '. There are two cases: 

1) If 7 = 7 ', we erase I and upgrade the weights of b and b' from ^fg to 


/—2a , \ ^ 2 a /c 

Vge + 2^ e Tk ’ 

k^l 

By (J3.27J) and the latter expression is bounded above in absolute value by 2e 2,3 . 

2) If 7 ^ 7', then we erase I and instead add two red links which connect between the endpoints 
of b and b'. Precisely, if b = (u,v) and b' = ( t,s ) where the both pairs {u,u} and {t, s} of dual 
vertices are recorded in the lexicographical order, then we add red links (u,t) and (v,s). In this 
way both red links lie on the dual lattice and have the same length k = 1,2,.... We associate the 
weight 77 to each of those links. 

Clearly after the above procedure is applied to all the excited intervals of T we end up with a 
connected edge self-avoiding polygon T which entirely lies on the dual lattice. In order to control 
the original weights u>(r) we over-count via ignoring the geometric constraints: from each vertex 
of the dual lattice one is permitted to grow up bonds in all 4 possible directions: either usual 
Ising bonds with weights e~ 2,3 or “red” bonds of lengths k = 1,2,... with the weights e 2ak rk 
respectively. Notice that any modified graph T contains at least 4 Ising bonds. Consequently, 


oer 


,2a |r| 


>(r)| 




Y 471 2e ~ 20 +Y ke 


2 a k 


Tk 


2=4 


(4.33) 


where the above sum is over the total number of bonds (and links ) of T. By (13.2611 and iTT 271) and 
in view of the possibility to control the smallness of 5n via (E23, we, given small q and a as in 
(THU . can always choose a large enough value of No, such that, 


2 e- 2 ?+ Yke 2ak r k 
1 


< 3e- 2 ^. 
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As a result the right hand side of (TOifl) is bounded above by 

(12e“ 2/3 ) 4 
1 - 12e~ 2 P' 

Recalling the notation e“ 2 ^ = y/g we, in view of (Id.201) . infer that the latter expression is much 
less than the value of a in l TP21 once q happens to be sufficiently small. 

In the sequel we shall assume that q ^ qo and N ^ No(q) are such that we actually have a 
strengthened version of 14.0511 : Set wo(r) = |w(T)|. Then, 

]T|r|e 2Q ! r l«;o(r) < a. (4.34) 

T90 

Indeed, (TO - ! follows by a straightforward adjustment of the arguments employed for the proof of 
Lemma E3 


4.4 Proof of Lemma 14.31 

Let q and No are fixed as in the proof of Lemma rOl and let {u>o(r) = |w(r)|} are the the absolute 
values of the weights of super-contours evaluated at such values of q and No(q). It is enough to 
check that there is q ^ q and No ^ No, such that for every N ^ No, the ( q,N ) super-contour 
weights {(w(T)} satisfy: 


Kr)|e^ er l r l ^ wo(r) and — 

z v 


yv 


(4.35) 


for every super-contour T and for each finite subset f CZ 2 . 

Of course, only the second inequality in (TT~35l) deserves to be checked. This is done by induction 
on the volume. Obviously, if the volume \V\ = 1, then we have ZylZy = 1. Suppose now that 
indeed 

^ < e\ dv \ , 

yh 

for all \V\ < K. We want to prove that this also holds when \V\ = K. In order to see that, observe 
that all the super-contours appearing in these two partition functions have interiors of volume at 
most K — 1. Introducing the sets and 5^ r _ 1 of all clusters made up of h-type, resp. v-type, 

super-contours having (total) interior of volume at most K — 1, and using the symmetry present 
in the model, we can write 


exp 




(£ 


Ecesb , \cnv\ 


xGV Z-iCgS k _ 1 
C3x 


~ 1 <j) T 


Z V 


exp 




\cnv\-i$T(C)) 

C3x 




Notice now that all the contours T appearing in the above partition functions have weights w(T) 
which, by the induction assumption and by the first of the inequalities in 14.3511 . satisfy: 


|ui(T)| ^ |w;(r)| e^'>' er 171 ^ wo(r). 

Therefore we can apply Lemma 14.11 Expanding the two partition functions and cancelling the 
terms involving clusters entirely contained inside V, we obtain the desired result, since by m 


vh 


v < 2a \dV\ 

.u O 


and 2a < 1 once, according to 632), q is not very close to 1. 
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5 Proofs of the main results 


In this section we complete the proofs of Theorem lL 2 l and Theorem Ol As in the proof of 
Lemma fl.ill we proceed to work within the range of parameters (q, N) which satisfy I4.H5I 


5.1 Contribution of long super-contours 

As before let V be a lattice torus of linear size L. Given a supercontour T £ define 


Gq t N,L(T) = 


J q,N,L 




F3F 


= i w{T)exp ^ -t> T (C) j . 


By (14.3511 and iTOoli . 

\G q , N ,Lp)\ <wo(r)e“l r l. 

Furthermore, by irm there exist constants C\ and C 2 such that 

Y. ^o(r)e“' r| < Cl L 2 e- afe ^«;o(r)e 2 a l r l 
|rj > k T30 

^ C2L 2 ae~ ak . 


(5.36) 


As a result, there exists C 3 < 00 , such that the contribution of super-contours T £ S k with 
|T| > CalogL to the partition function Z q ^,L is, uniformly in L , negligible. The same argument 
applies, of course, in the case of vertical super-contours T £ 5|. 


5.2 Proof of Theorem 11.21 

The first statement of Theorem 11.21 follows immediately from results of Gruber and Kunz 0. 

Assume now that the parameters (q, N ) satisfy (14.351) . By 14.3411 we may exclude long winding 
contours. Thus, the only thing remaining to be done in order to complete the proof of Theorem ll.2l 
is to estimate the probability that a given site, say 0 , belongs to the interior of some short non¬ 
winding contour. In view of Lemma 14.11 the probability that 0 belongs to the interior of such a 
super-contour can then be written as 

= ^ a(r)exp ( £S T (C)) . (5.37) 

roo ^q,N,v roo y c^r J 

where TOO means that 0 is in the interior of the super-contour T. By the latter expression 

is bounded above by 

YMr)e alrl < ^|rK(r)e a|r| , 

roo rso 

The claim of the Theorem follows now from (irm 


5.3 Infinite volume states 

Let Aoo be the set of all such coverings u> of 1? by horizontal and vertical rods (we colour monomers 
as well), which contain only finite contours. Of course, for every ui the notion of the exteriour colour 
x(u;) = h or u is well defined. By a straightforward application of Lemma 14.11 
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Theorem 5.1 There exists qo > 0 such that for every q ^ qo one can find No = -/Vo( 5 ) which 
enjoys the following property: For every N ^ No there exists a unique infinite volume Gibbs state 
h’qN (respectively p ? qlv / such that 

(A 00 ; X (Z) = h) = 1 

(respectively N (AooWfix) = v) = 1) . Furthermore, let 7 be the (random,) set of of all the 
exteriour contours of ui and, given a finite domain A C Z 2 , let 7 A = ( 71 ,..., 7 n ) be a fixed 
compatible set of exteriour contours, such that each 7 *, intersects A, A fl 7 *, 0. Then, 


(l A c 7) = w (Ea) ex P 

- £ <h T (C)| 

I A ~7 \ 

V C^La J 


(5.38) 


where the above sum is over all compatible collections T A = (IA,... , r m ) of super-contours satis¬ 
fying: 

V l = 1 ,... , m3 k such that 7 ^ £ T; and U 7 fcCur;. 

Formulas and (TOol) readily imply that ia q N has an exponential clustering property: Given 

two disjoint boxes Ai and A 2 and two fixed compatible collections 7 a andy A with 7 A C Aj; k = 
1 , 2 , the following bound holds: 

log N ( 7Al c l'^ A 2 C l) 

Tq, N (7A! C 7 ) /^iv (7 a 2 C 7) 

where d(Ai, A 2 ) is a distance (say h) and ci and C 2 are two positive constants which depend only 
on q and N. 


< c 1 |A 1 ||A 2 |e- aC2d(Al ’ A2) , (5.39) 


5.4 Boundary surface tension 

Consider vertical and horizontal intervals 

r k = (1/2,1/2) + {(0, 0), (0, 1 ),..., ( 0 , k - 1)} and j{) = (1/2,1/2) + {(0, 0 ),... (k - 1 , 0)}. 

By construction, both Jf and Jf are linear segments on the dual lattice (1/2,1/2) + Z 2 . Given a 
rod / = (ui,..., u n ) C Z 2 let us say that / intersects Jf; I fl Jf y^ 0 if 

JJnint(Uj =1 B 1 (« fc ))#0, 

where Bi{u) = u+ [—1/2,1/2] x [—1/2,1/2] C R 2 and for a bounded set let 2 the symbol int(A) 
stands for its R 2 -interiour. In a similar fashion we define lT\J k y^ 0. Notice that monomers cannot 
intersect Jf or J k . Also, with such a definition, J k cannot be intersected by a vertical rod and, 
accordingly, J k cannot be intersected by a horizontal one. 

Given a Z 2 tiling ui £ Aoo let us say that the event {Jf fl ui = 0} ( respectively { J k fl ui = 0}) 
occurs if Jf (respectively J k ) does not intersect any of the rods of u>. 

We define two types of boundary surface tensions: 

r q ,N = - lim j log fi qN (Jfc n <J = 0) = - lira j log/^ jJV (J k fl w = 0) , (5.40) 

#C ^ OO (V (C ^ OO /v 

and 

(q,N = - lim \ log ii h qN (j/nu = 0) = - lim j log ( J k D u = 0), (5.41) 

In both cases the fact that the corresponding quantities are well defined follows from standard sub¬ 
additivity arguments based on the exponential clustering property and on the 7r/2-rotational 

symmetry between the vertical and horizontal states. 
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Lemma 5.2 For any q sufficiently small there exists No = No(q) such that for every N ^ Nq, 


T q,N > %q,N- 


(5.42) 


Proof. Since we do not try to prove the lemma in the whole range of entropy driven symmetry 
breaking, the poof boils down to a crude perturbative argument. We start with a lower bound on 
r q^N ■ Fix 71 ,..., to be the set of all exteriour contours of lo wI nch intersect • 11 y n> and 

lEHTl) . 


Mg,TV FI Jk 



< (2q) |J ^ U7i| . 


(5.43) 


It remains, therefore, to derive an upper bound on 


^,jv(l^\urnl<fc/2). 

By a straightforward modification of the over-counting argument employed in the proof of Lemma l4.2l 
we infer from 115.3811 that for a given collection 71 ,... ,7 n of exteriour contours, 


7i»--->7n) < exp (- 2 / 3 ^diam( 7 j) 


1 


where, as before, e 2,3 = ffig and g is related to q via 13.2011 . Elementary combinatorics leads then 
to the following conclusion: If q is sufficiently small and N ^ No(q), then 


/£*(!■# \U7l|<fc/2) 



(5.44) 


Combining 15.4311 and 15.4411 we arrive to the following lower bound on t^jv: 

Tq ,n > - ^ logg. 


(5.45) 


In order to derive a complementary upper bound on jv notice that on the level of events (under 
the vertical state ti v q N ), 

{cu fl = 0} 3 {V 7 exteriour contour of w Jf D int( 7 ) = 0} . 


Indeed, by the definition Jf can be intersected only by horizontal rods. Let us say that a super- 

i 

contour F is intersection incompatible with J£; F 7 ^ Jf, if T contains a contour 7 , such that 
Jf fl int( 7 ) ^ 0. Then, by Lemma E3~T1 


/ 


Mg.JV (w n Jfc = 0 ) ^ exp 


\ 


—ak 


- Y, ! $T ( C )I > e 

V c/jj ] 

Consequently, ^ a and, in view of 15.4511 and (IPH . the proof of Lemma O is concluded. 


5.5 Sketch of a proof of Theorem 11.31 

Consider boxes Vn with periodic boundary conditions. As before we continue to ignore winding 
super-contours. In particular the notion of exteriour colour is always well defined. Let, therefore, 
fc 6r an d -^n’fe 01 be the partition functions of rod tilings of Vn with the exteriour colour being 
fixed as h (respectively v). By (I5.3(il) and Lemma rm 


log 


r^h, per 

^n,k 


rrV, per 

^n,k 


^ c%n 2 e 


—c^an 


(5.46) 
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Finally, let and p^ p k el be the corresponding Gibbs states. 

The partition functions Z^’ k and of the (exteriour colour) horizontal and vertical tilings 
of Vjr with free boundary conditions are related to Z^’ pel and as follows: Set 

J n,k = {* = (*i»*2) € V£ : ii = 0} and J^ k = {i = {i u i 2 ) S V£ : i 2 = 0} . 


Then, 


and, respectively, 


yh,f 
^n,k 
yh, per 
Z n,fe 


lyV, f 

yv,per 

^n,k 


VnT ( u n Jn,k = 0 ; w n = 0), 
<r( 5nJ 4 = 0;SnJ^ = 0). 


By (15.4011 and 15.4111 the latter probabilities are logarithmically asymptotic to 


exp (—n((2fc 2 + l)r g ,jv + (2*4 + l)£ 9 ,iv)) and exp (-n((2fci + 1 )t 9i jv + (2fc 2 + l)£g,Jv)) 
respectively. The claim of Theorem o follows now from and Lemma IQ 
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